Lifetime measurements of the 5d states of rubidium 
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We present lifetime measurements of the 5-D3/2 and 5-D5/2 states of rubidium using the time 
correlated single photon counting method. We perform the experiment in a magneto-optical trap of 
*^Rb atoms using a two-step excitation with the trap laser at 780 nm as the first step. We record 
the 761.9 nm fluorescence from the decay of the 5-D3/2 state to the 5Pi/2 state, and measure the 
lifetime of the 5-D3/2 state r = 246.3(1.6) ns. We record the 420.2 nm fluorescence from the cascade 
decay of the 5-D5/2 state to the 5Si/2 state through the 6P3/2 state, and extract the lifetime of the 
5D5/2 state r = 238.5(2.3) ns. 



I. INTRODUCTION 

Lifetime measurements are fundamental for the un- 
derstanding of atomic structure. The precisely mea- 
sured quantities of transition energies and excited state 
lifetimes are important tests of atomic structure theory 
through comparisons with the values derived from the 
calculated energies and transition matrix elements 
Alkali atoms and alkali-like ions are the benchmark sys- 
tems where measurements and calculations have reached 
the highest development because the simple single-core 
electron structure facilitates the application of advanced 
calculation methods and current simple laser cooling 
techniques. 

A large number of lifetime measurements exist in the 
alkali atoms and alkali-like ions for the s and p levels, 
but there are not as many high precision measurements 
on the d and higher angular momentum states. These 
states are becoming important not only in the study of 
fundamental symmetries [1, H, 0, Q , but also in quantum 
information science as they are used for qubit manipula- 
tion in ion traps 0] ■ 

The d states of rubidium are interesting from the calcu- 
lation point of view as they remain in the domain of non- 
relativistic physics, but they are thoroughly affected by 
correlation corrections. Safronova et al. 8] have shown 
the important role of high order corrections, up to third 
order, in calculations that use many body perturbation 
theory (MBPT) of lifetimes of these states. We present 
in this paper a precise measurement of the and 
5D5/2 states of rubidium. Our group has studied the 
equivalent states in francium, the heaviest alkali atom, 
and achieved precisions of 0.4% and 4.3% for the 
and 71)5/2 states, respectively @. 

Previous experimental work on the lifetime of the Rb 
5D3/2 state and the whole 5c? manifold [llj achieved 
a precision long surpassed by atomic calculations. This 
paper shows an improvement by more than a factor of 
twenty in the precision of the 51)3/2 state lifetime mea- 
surement and a precise determination of the bDc^ state 
lifetime. These improvements will hopefully trigger im- 
proved calculations as the precision of this work (^ 1%) 
is better than the current estimates in the theory (j^ 5%) 



This paper has five sections. Sec. [IT] presents the the- 
oretical background. Sec. IIIII explains the experimental 
method and apparatus. Sec. IIVI and fVl presents the data 
analysis for the measurements of the 5D3/2 and 51)5/2 
state lifetimes, respectively. Sec. IVII gives a conclusion. 



II. THEORETICAL BACKGROUND 

Suppose an atom is in an excited state |e), and it de- 
cays to a set of lower states the lifetime of this 
excited state, Tg, is Q 
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where c is the speed of light, hio^i- is the transition energy 
between |e) and a is the fine structure constant, J is 
the angular momentum associated with each state, and 
(Je||r||J/.) is the reduced transition matrix element. 

Equation ([T]) shows that lifetime measurements probe 
the overlap of wavefunctions of electron states with dif- 
ferent quantum numbers at large distance. This is in 
contrast with precise energy measurements which probe 
the overall distribution of electron wavefunction, or the 
hyperfine splittings that are sensitive to the short range 
of the electronic wavefunctions. 

Theoretical calculations aim to first find out the correct 
set of eigenfunctions of the atomic system Hamiltonian, 
and then calculate the physical properties using those 
results. In relativistic many-body calculations, there are 
different methods to obtain the eigenfunctions of the sys- 
tem. The most straightforward way is to start with a 
Dirac-Hartree-Fock (DHF) wavefunction dV''"')) and it- 
erate the perturbation theory order by order. This is the 
relativistic many-body perturbation theory. This method 
usually stops at the third order due to the large calcula- 
tion load. This limits the precision of the results that can 
also suffer when the mean field ground state is not a good 
starting point. An improvement is to use a wavefunction 
including some excitations to the DHF wavefunction [l|: 

|V) = (l + ^l+52+53 + ...)|V^(*"), (2) 

where j')/''^''') is the DHF wavefunction and Si is the i- 
body excitation operator. The method involving only S\ 
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and 5*2 operators in Eq. ([2]) is the singles-doubles (SD) 
method. The so-called all-order method includes 5*1, S2 
and the dominant parts of the S'3 (see Ref. [l^ for an 
extended explanation of the method). 

The calculation of transition matrix elements limits 
the precise prediction of the lifetimes. For the all-order 
method, this problem becomes serious for some levels 
with large correlation effects (such as d levels), where 
the wavefunction in the form of Eq. ([2]) is a bad starting 
point of calculations. Theorists use a semi-empirical scal- 
ing scheme for certain wavefunction coefficients to solve 
the problem. For the rubidium 5d states, the lifetimes 
calculated with and without scaling have a relative dif- 
ference of 20%, but it is hard to assert the correctness of 
the scaling due to high uncertainties of past experimental 
results [8i]. 

Lifetime measurements help on the extraction of tran- 
sition matrix elements. If there is only one decay 
channel, we can extract the absolute value of the re- 
duced transition matrix element directly from the life- 
time data. For example, by measuring the lifetimes of 
the excited states, Ref. [13] extracts the transition ma- 
trix elements i(7Pi/2||r||75i/2)| and |(7P3/2||r||7^i/2)| of 
Fr, and Ref. [ij extracts the transition matrix element 
|(5£'5/2||?'||6-P3/2)| of Cs. The case of a multi-channel de- 
cay sometimes requires a combination of various exper- 
imental approaches to extract the transition matrix ele- 
ments [l| , including the measurement of static scalar po- 
larizability TB'l, and nonrcsonant two-photon, two-color 
linear depolarization spectrum ,.16] . 



III. EXPERIMENTAL METHOD AND SETUP 

We adapt the time correlated single photon counting 
method in this experiment with a cold sample of *^Rb 
atoms (temperature less than 100 /xK) in a magneto- 
optical trap (MOT). After exciting the atoms to the bd 
states, we turn off the excitation beams and record the 
delay between the fluorescence and a fixed time reference. 
Here, the single photon counting technique refers to the 
fact that we record at most one photon in one duty cy- 
cle. It is possible for the detector to detect more than 
one photon in one cycle, but once the electronics record 
one signal, they will not accept another until a new cycle 
begins. This method works best when the probability 
of detecting more than one photon in one cycle is very 
small, which in turn keeps the corrections low (see the 
pulse pileup corrections in the systematic study of Sec. 
lIVp . In our experiments, we detect one photon in about 
one hundred cycles. 

We use a rubidium dispenser as the atomic vapor 
source and the MOT resides inside a 15 cm radius spher- 
ical chamber with the vacumm pressure of 10~^° Torr. 
A pair of anti-Helmholtz coils provides a magnetic gradi- 
ent of 6 G/cm and three pairs of Helmholtz coils provide 
the fine tuning of the magnetic environment. A Coher- 
ent 899-01 Ti:Sapphire laser with linewidth better than 




FIG. 1: Rb energy levels involved in the experiment. 



100 kHz provides three pairs of MOT trapping beams 
with intensity of 8 mW/cm^, and the laser is red de- 
tuned from the 55*1/2, F = 2 ^ 5P3/2,F — 3 transition 
by approximately 20 MHz. A Toptica SCllO laser pro- 
vides the repumper beam with intensity of 3 mW/cm'^ 
and it is on resonance with the transition 551/2, F = I 
5P3/2, F = 2. We capture about 10^ atoms in the MOT 
with diameter of 600 /im and peak density of around 10^ 
cm~'^. We use two charge-coupled device (CCD) cameras 
to monitor the fluorescence of the MOT in two perpen- 
dicular directions. 

We list the relevant energy levels of *^Rb for this ex- 
periment in Fig. [1] We use a two-step transition to reach 
the 5d states, where the trapping beam of the MOT en- 
ables the first step and the 5P3/2,F = 3 state is the 
intermediate state. A SDL TC40-D laser with linewidth 
of 5 MHz provides the probe beam to reach the 5d states. 
We inject the probe beam to the MOT region through a 
single mode fiber, which sets the waist (1/e^ power) to 
1.2 mm. The power of the probe beam is 0.5 mW for the 
excitation to the 5D^/2, F = 4 state, and 1.0 mW for the 
excitation to the 5I?3/2, F — 3 state. 

We lock the frequency of the trapping beam using the 
Pound-Drever-Hall method with saturation spectroscopy 
in a Rb cell. We send part of the frequency modulated 
light employed on this lock to an independent rubidium 
glass cell, where this light overlaps with a small laser 
beam taken from the probe beam. We monitor the ab- 
sorption of the 780 nm light as a function of the frequency 
of the 776 nm laser, and the intermodulation of the side- 
bands yields error-signal like features that we use to lock 
the frequency of the probe beam on resonance [l^, [l^ . 

We use a cycle of 10 ^s length, and have two different 
schemes for photon detection and time control for the 
two different measurements. In the 5I?3/2 state lifetime 
measurement, we put a 760 nm interference filter with 
bandwidth of 10 nm (Andover 760FS10-25) in front of 
the detector, a Hamamatsu R636 photomultiplier tube 
(PMT) with quantum efficiency of 10% at this wave- 
length. Since a large amount of 780 nm photons from 
the scattered trapping beam pass through the filter, we 
turn the trapping beam off after the excitation phase 
to decrease the background. We use two acousto-optical 
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modulators (AOM) to turn on and off the trapping beam 
and the probe beam. In the 5D5/2 state hfetime measure- 
ment, we use a 420 nm interference filter with a band- 
width of 10 nm (Andover 420FS10-25) in front of the de- 
tector, a PerkinElmer C962 channel PMT with quantum 
efficiency of 15% at this wavelength. The trapping beam 
has a negligible contribution to the background in this 
case. We use a Gsanger LM0202P electro-optical modu- 
lator (EOM) and an AOM to chop the probe beam for 
this experiment. The turn off ratio of EOM and AOM is 
better than 30:1 in 30 ns. Fig.[2]shows the time sequence 
of the pulses. 
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FIG. 3: The setup of the detection and signal process. 
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FIG. 2: Time sequences in the hfetime measurement. Each 
cycle has a period of 10 ^s. 



The PMT is about 35 cm away from the MOT and 
we use a lOX Computar Macro Zoom lens in front of it 
for light collection. Two synchronized Stanford Research 
DG535 pulse generators, which have a 5 ps delay resolu- 
tion and 50 ps rms jitter, provide all the time references 
in the signal process. The signal from the PMT goes 
through the amplifiers first. In the 5Z?3/2 {50^/2) life- 
time measurement, we amplify the signal from the Hama- 
mastsu R636 PMT (channel PMT) 64 (4) times using 
an EG&G AN106/n plus an ANlOl/n DC amplifier (an 
AN106/n amplifier alone). The output goes through an 
Ortec 583 constant fraction discriminator (CFD) and a 
Lecroy 7126 level translator. The level translator con- 
verts the input signal to ECL, TTL and NIM outputs. 
We send the output of the NIM signal to a Stanford Re- 
search SR430 multi-channel scaler to monitor the photon 
counting histogram during the experiment. We send the 
ECL signal as a start pulse to a Lecroy 3377 time-to- 
digital converter (TDC), which has a resolution of 0.5 ns 
and is triggered by the falling edge of the input pulse. 
The TDC measures the delay between the observed pho- 
ton and the fixed pulse given by the pulse generator (see 
Fig. [2] and Fig. [S]). The output of the TDC goes to a 
Lecroy 4302 memory and we read out the results through 
a Lecroy 8901A GPIB interface. Fig.[3]shows the diagram 
of the signal process. 



IV. LIFETIME OF THE 5D3/2 STATE 

We excite the atoms to the 50^/2 state and about 50% 
of the atoms decay into the 5Pi/2 state with fluorescence 
at 761.9 nm, where we get the branching ratios using the 
transition matrix elements from Ref. [8|. We record the 
fluorescence and accumulate data until the peak count is 
more than 1000 for each data set. This usually takes half 
an hour to 45 minutes. The ratio of peak signal to back- 
ground is generally around 40, and we take an additional 
data set of background for roughly the same time. The 
upper plot of Fig. 2] shows a typical data set after sub- 
straction of background. There is another possible way 
to observe the decay through the detection of 420 nm flu- 
orescence that comes from the 6p manifold (see Fig. [T]). 
This has a great background advantage, but requires flt- 
ting to a sum of three exponentials, which challenges our 
fltting process (see explanations in the Sec. |V|. The re- 
sults that we present here are consistent with a prelimi- 
nary analysis of the blue decay. 

A. Decay model 

The data in Fig.|3]shows an increase of the fluorescence 
due to the excitation of the system by the probe beam 
and the decay after turning off the probe beam. We 
qualitatively understand the rising part by treating this 
many-level system as a four-level system. We denote the 
ground state, the 5P3/2 state and the 5D3/2 state as levels 
1, 2 and 3, and treat all other states that the 50^/2 state 
decays to as level 4. The excitation rate from level i to 
j is Aij , the spontaneous decay rate of level k is Rk and 
the branching ratio of decay channel from fc to Hs 'jki- 
Because the system is already in the steady state before 
we turn on the weak probe beam, we take the second step 
transition as a small perturbation to the system and get 
the rate equations: 

N3 = A23N2 - R3N3 

N4 = I34R3N3 ~ RiNi. (3) 

The initial condition is A^i(O) = N2{0) = No/2, where 
A'o is the total number of atoms and the time origin is the 
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FIG. 4: Upper plot: time distribution of 761.9 nm photons 
(black points) in the measurement of the 5D3/2 state with 
the fitting function (white line). Lower plot: normalized fit 
residuals. We have subtracted the background in the decay 
curve and we use the statistical error of the data to normalize 
the fit residuals. 

time we turn on the probe beam. This gives the number 
of atoms in level 3 as 

^^3W = ^^[l-exp(-i?3i)],t<ti, (4) 
its ^ 

where ti is the time we turn off the probe beam. This 
qualitatively agrees with the rising behavior of the signal. 

The number of atoms in this level begins to decay when 
we turn off the probe beam. For a perfect turnoff at time 
ti, we get the decaying behavior: 

Nait) = Nsih) exp [-R^it ~ h)],t > h. (5) 

The intensity of the fluorescence from the 5I?3/2 state 
to the 5Pi/2 state is proportional to the atom number in 
the 5D3/2 state and Eq. ([5]) shows that the decaying part 
of the signal is independent of the excitation part except 
for a coefficient. 

B. Fitting process 

We use the Levenberg-Marquardt (L-M) algorithm [20j 
and flt the data to the function 

y = Aexp{^t/T) + B, (6) 



where the r is the 5D3/2 state lifetime and B describes 
the background. The fitting is a procedure to minimize 

X^^Y.''^{y.-Vf)\ (7) 

i 

where yi and y/ are the experimental and fitting value for 
the zth data point, and Wi is the inverse of the square of 
the statistical error of the data assuming the data follow 
the Poisson distribution. 

We use two criteria on the choices of the range of the 
data to fit. First, we choose the starting point of the 
fitting range away from the turning off point to avoid 
pulse turnoff effects. Second, we need to make sure that 
the fitting results are stable when we vary the starting 
and end points of the fitting range. This fitting process 
gives us an average of reduced — 1-03 for all the data 
sets we have taken. The lower plot of Fig. 2] shows the 
normalized fit residuals of the upper plot data. 

C. Statistical error 

We have taken 28 sets of data with the same experi- 
mental conditions, with each individual data set fitting 
error, ct^, of around 1%. Considering this sample of 
28 data points, we obtain a weighted mean lifetime of 
246.5 ns, with the error given by (X) 1/"'?)^^''^: which is 
0.3 ns. We choose to scale the uncertainty by a factor 

m 

S^[-X^/{N-l)]\ (8) 

where is defined in Eq. ([7]), and N is the number of 
data points. In this way, the mean value does not change, 
but it helps to solve the problem of underestimation of 
the error on some data sets. We get S — 2.1, and the 
statistical error as 0.6 ns. 



D. Systematic effects 

We focus on the following systematic effects: 
Time calibration and height nonuniformity of TDC. 
We send two pulses with a fixed delay from the DG535 
pulse generator as the start and end pulses to the TDC. 
By comparing the readout of the TDC with the calibrated 
delay, we get a nonlinearity of the TDC time calibra- 
tion of less than 0.01%, which in turn gives an error on 
the lifetime measurement of less than 0.01%. To check 
the height uniformity of the TDC channels, we change 
the source of the photon to the scattering of laser light 
on a paper. We control the power of the laser so that 
the recording rate is similar to that in the lifetime ex- 
periments, and accumulate data until the count of each 
channel reaches 100. We fit the results with a linear func- 
tion (fits to higher order polynomials yield the same con- 
clusion) and get a slope of 0.1 count per 1000 channels. 
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which gives an error of 0.01% on the hfetime measure- 
ment. 

Pulse pileup correction. There is a difference between 
the number of the signals detected and signals recorded 
in the single photon counting method, this leads to pulse 
pileup correction. Suppose we have Ne cycles of excita- 
tion, and get Nj counts on the jth channel of the TDC, 
then the pulse pileup corrected result of the ith channel 



N' 



1-^ 



(9) 



This gives us a correction of around -0.1% for our mea- 
surements. 

Quantum heats. Quantum beats come from the inter- 
ference of the decay paths from several coherently excited 
states to the same lower state. By considering this effect, 
we rewrite the Eq. © as [13] 



y = A exp(-t/T) [1 + ^ C, cos(cJii)] 



(10) 



where the sum is over all possible interference paths and 
hbJi = AEi is the energy difference of two excited states 
involved in the process. This has been an important issue 
when the probe beam has a large bandwidth [23|. 

We first consider the possibility of the quantum beats 
from the hyperfine structure. The energy difference be- 
tween 5Z?3/2, F = 2 and F = 3 sates is around 40 MHz 
[To| . This separation is much larger than the probe laser 
linewidth. The natural linewidth of the 5D3/2 state is 
0.6 MHz, considering the power broadening effect, the 
transition linewidth is at most 6 MHz. A sharp turnoff 
process could generate a large bandwidth in frequency 
domain. We decompose the probe pulse in frequency 
range by the fast Fourier transform (FFT) , and find that 
the high frequency components (higher than 40 MHz) 
only take about 1/200 of the total power. In conclusion, 
the probability to coherently excite the hyperfine states 
is negligible. The FFT of the fitting residuals does not 
show any peaks around the frequency of the hyperfine 
splitting either. 

We position the atoms at the zero of the magnetic field 
in the MOT, but the finite extent of the cloud sample re- 
gions with a magnetic field. The Zeeman splitting intro- 
duced by this magnetic field could be a source of quan- 
tum beats, because we always leave the gradient magnetic 
field for the MOT on during the experiment. However, 
the overlapped a+ and a— trapping beams tend to av- 
erage out this effect. 

We estimate an upper bound of the effect due to the 
Zeeman splitting following the treatment in Ref . . We 
have huj = AmpgpiiBB, where Amp ~ 2, gp ^ 0.4 for 
5D^/2,F = 3 state, and /is is the Bohr magneton. The 
size of the MOT is around 600 fim and the magnetic gra- 
dient is around 6 G/cm, we use the maximum magnetic 
field 0.18 G in that region for the calculation. This gives 



us a maximum beat frequency of 202 kHz, which sets an 
upper bound of this error of 0.15%. 

We also look for possible effects of the magnetic field 
by changing the magnetic field gradient from 5 G/cm to 
10 G/cm and do not find any systematic changes. Fig. [5] 
shows the results, where the smaller error bars in the 
points with 4.9 G/cm and 9.8 G/cm are due to more 
data sets under those conditions in this systematic study. 
Studies of linear correlation coefficient of the data points 
show that they are highly uncorrelated. We get an error 
of 0.16% by this effect, which is the standard deviation of 
the data points taken in this study. We use this measured 
bound of 0.16% for possible magnetic effects that include 
the quantum beats. 
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FIG. 5; The systematic study of the effect of the magnetic 
field gradient on the 51)3/2 lifetime measurement. 

Radiation trapping. Radiation trapping comes from 
the reabsorption of the emitted photon by the sample 
itself. This effect depends on the length and density of 
the sample and could increase the measured lifetime sub- 
stantially when there are many atoms [131 • The MOT has 
the advantage of its size, which makes this effect much 
smaller compared with the experiment in the cell. De- 
note I and n as the length and density of the sample, and 
a as the atomic absorption coefficient. Ref. [1^ gives an 
expression for this effect: 



where C = Ian and Atan(A) = C. 
length, Eq. pT|) becomes 



(11) 



In the limit of small 



tC — rlan, 



(12) 



where the change of lifetime depends linearly on density. 
We have changed the density of our sample by a factor 
of six to look for this effect, and do not observe any sys- 
tematic changes. The error due to this effect is smaller 
than 0.1%. 

Other systematics. The pulse length and power of the 
probe beam are the parameters controlling the number of 
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TABLE I: Error budget of the 5-D3/2 state lifetime measure- 
ment. 



Source 



Correction (%) Error(%) 



Statistical 
Time calibration 
TDC nonuniformity 
Pulse pileup 

Quantum beats and magnetic field 
Radiation trapping 
Other Systematics 



-0.1 



±0.25 

< ±0.01 
±0.01 

±0.16 

< ±0.10 
±0.60 



Total 


±0.66 


TABLE IL Comparison of the measured lifetime of the 5D3/2 


state with the previous work. 






Experiment This work 


246.3±1.6 


Tai et al. ffO] 


205±40 


Theory Theodosiou [26] 


240 


Safronova et al. [8| 


243 



atoms in the excited states. By changing the pulse length 
from 320 ns to 720 ns, and also the power from 0.5 mW 
to 2.4 mW, we do not observe systematic changes. We 
set the error on those effects as 0.6%, which we obtain 
in the same sprit as in the study of the magnetic field 
effects. 




400 600 800 

time (ns) 

FIG. 6: Upper plot: time distribution (black points) of the 
blue photon in the measurement of the 5D5/2 state with 
the fitting function (white line). Lower plot: normalized fit 
residues, where the normalization factor is the statistical error 
of the data. 



E. Summary and comparison 



LIFETIME OF THE 5D5/2 STATE 



We list the error budget of the 5D3 /2 lifetime measure- 
ment in table |T1 and compare our result with the previous 
experimental and theoretical work in table [TTl On the ex- 
perimental part, Ref. 10] made use of the Hanle effect 
to extract the lifetime. In that case, they could not op- 
tically excite the atoms to the 5D3/2 state, instead, they 
optically excited the atoms to a higher p state and popu- 
lated the desired level by spontaneous emission. On the 
theoretical part, Ref. [20| used the quantum defect the- 
ory with a realistic potential, and Ref. [1] used the scaled 
all-order method. 

The result of Ref. [1] sits on the edge of 2a range of 
our data and that of Ref. [2^ sits on the edge of the 4cr 
range. For more detailed comparison, we need to know 
the range of validity of the theoretical calculation. A 
possible way to define the error of the all-order method 
calculation is comparing its results with third order many 
body perturbation theory. In this way, the error of the 
all-order method calculation for the Rb bd states lifetime 



is around 5% |27[, and our result of the 5D3/2 state life- 
time stays within its prediction. Our result Icr range has 
a small overlap with the Icr range of the result of Ref. p^ . 
which has a much bigger error bar. 



About 30% of the atoms in the 5D5/2 state decay to 
the 6P3/2 state with fluorescence at 5 /xm, and about 
25% of the atoms in the 6P3/2 state decay to the ground 
state with fluorescence at 420.2 nm, where we get the 
branching ratios in the same way as in Sec. lIVI We record 
this blue fluorescence and accumulate data until the peak 
count is more than 1000 for each data set. The upper plot 
of Fig. [6] shows a typical data set. Since we detect the 
blue photon, the background is very small (less than 1 
count per channel in half an hour), and we do not have 
to take an additional background data. 

Most of the data process and systematic effects studies 
in this measurement are similar to those in the measure- 
ment of the 51)3/2 state, however, there are some new 
features in the analysis. 



A. Decay model 

We use the same model in the previous section. The 
differences are that we change the 5Z?3/2 state to the 
51)5/2 state and specify level 4 as the 6P3/2 state, be- 
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cause the fluorescence we observe in this experiment is 
proportional to number of atoms in the 6P3/2 state. Us- 
ing the same notation as in Sec. IIV A[ we get the rising 
and decaying part of the signal in this experiment as 



i?4 — -R3 2 

1 — exp(— i?3t) 



1 — exp(— i?4t) 



i?3 



-R4 



. t <ti 



N4h) - 
734^3 

i?4 - R; 



J34R3 
Ra — R; 



exp [— i?4(t — ti)] 



-Naih) exp [-R3{t - ti)],t > ii(13) 



This shows that the curve observed in this experiment 



5/2 



includes the information of the lifetimes of the 5D, 
state and the 6P3/2 state. 



B. Fitting process 

The ratio between coefficients of the two exponentials 
depends on the lifetimes and initial atom number at the 
turnoff point of both states (see Eq. (|13p V If we want to 
fit the decay curve according to the function in Eq. (fT3|) . 
we have to evaluate the atom number by doing the con- 
volution of the excitation pulse with the formula that 
determines population in the 6P3/2 state. A simple way 
to avoid this process is to leave both coefficients as free 
parameters and fit the decay signal away from the turnoff 
point The existence of a stable range of fitting re- 
sults when we change the starting and end points of the 
fitting range is a test of the validity of this treatment. 

We use the following fitting function to describe the 
data 



y = Ai exp(-t/ ) -t- A2 exp(-i/ tqp^ 



B, (14) 



where t^d.^/^ ^'^d tqp^^^ are the lifetimes of the 5D5/2 and 
6P3/2 states, respectively and B is a possible background. 

If we leave both lifetimes as free parameters, the L-M 
algorithm does not give reliable results (it tends to equate 
both lifetimes). Therefore, we fix the value of TePg^^ in 
the fitting process to solve this problem, and the value 
we use is the weighted average t^^p^^^— 112. ^±.1.1 ns from 
previous experimental work listed in Ref. 11]. Our fitting 
process gives an averaged reduced of 1.09 for all the 
data sets we have taken. The lower plot of Fig. [6] shows 
the normalized fit residuals by comparing the raw data 
and fit result in the upper plot. 



C. Statistical error 

We have taken 10 sets of data with the same experi- 
mental conditions, with each individual data set fitting 
error, at, of around 1%. Following the same procedure in 
the Sec. IIVI we get a weighted mean lifetime of 238.6 ns, 



error of the mean as 0.8 ns, and the scaling factor 5 = 1.2. 
We choose the scaled error of the mean as the statistical 
error, which is 1.0 ns. 



D. Systematic effects 

We have studied all systematic effects listed in the mea- 
surement of the 51)3/2 state, which shows similar results. 
We focus on the following new situations in this measure- 
ment: 

Quantum beats. It is possible to coherently excite the 
atoms to the Zeeman sublevels of the 5D5 /2 state, but to 
observe the quantum beats signal for this cascade decay, 
we also need to prepare the coherence on the 6P3/2 state, 
for example, by means of selecting polarization on the 
fluorescence [28|. In our experiment, we do not attempt 
to observe the coherence in the intermediate level, and do 
not expect observable quantum beats effects. We check 
the Fourier transform of the fitting residuals and do not 
find peaks. We also check the magnetic field effects on 
the lifetime measurement and do not see any systematic 
effects. We set the error as 0.44%, which is defined in 
same way of that in the systematic study in Sec. IIVI 

Bayesian error. In the fitting process, we do not treat 
and T5£)j.^2 as independent parameters because we 
fix the T6P3^2 by its experimental value. Hence there is 
correlation between the values of tqp^^^ and 
the fitting process. We express the probability of finding 



T5D5/2 equal to t{ as 



(15) 



where denotes the possible value of tqp^^^ and -P(t{ jrj) 
is the conditional probability of getting tsd^^^ = t[ when 
'''6P3/2 — '^2- Eq. (|15p shows that the error of Tgpg^j value 
propagates to that of T^Dr^j-i value through the correla- 
tion, and this is the Bayesian error. 

Following the treatment of Ref. [29|, we change the 
value of Tgpgyj from 108 ns to 118 ns with an increase of 
0.5 ns each step, and do the fitting process to extract a 
new T5£)5^2 for each new value of TePg^^. This scanning 
range has covered 3cr range of the experimental result of 
T6P3/2 • We find that a linear function describes well the 
dependance of the fitting result of '^'^ value of 



3/2 



n(T2) = Tl -f a(T2 - T2), 



(16) 



where the notation is the same as in Eq. ([15)1 . Hence we 
get the Bayesian error as \a\arep^^^. For our data sets, 
the analysis gives a as -0.3, and the Bayesian error as 
0.5 ns. 

Other systematics. We study the pulse length and 
power effects, and do not find any systematic changes 
in the lifetime values. The limit set on this error due to 
these effects is 0.73%. 
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TABLE III: Error budget of the 5-D5/2 state lifetime measure- 
ment. 



Source 



Correction (%) 



Error (%) 



Statistical 
Time calibration 
TDC nonuniformity 
Pulse pileup 
Magnetic filed 
Radiation trapping 
Bayesian 

Other Systematics 
Total 



-0.05 



±0.40 

< ±0.01 
±0.01 

±0.44 

< ±0.1 
±0.21 
±0.73 

±0.96 



TABLE IV: Comparison of the measured lifetime of the 5-D5/2 
state with the previous work. 



(ns) 



Experiment This work 238.5±2.3 

Marek and Munster " [llj 230±23 

Theory Theodosiou [26] 231 

Safronova et al. \8} 235 

"This measurement gave a value for the whole 5d manifold. 



E. Summary and Comparison 

We list the error budget of the 5D^/2 state lifetime 
measurement in the table[lll]and compare our result with 
the previous experimental and theoretical work in table 
rvl The result of Ref. B sits within the 2a range of our 
data, and that of Ref. [2g| sits close to the edge of the 3a 
range. Since the error of the all-order method estimation 
on this level is 5% (see the Sec. lIVE|) . the result of Ref. [i| 
is consistent with our result. On the experimental part, 
as far as we know, there has been no explicit experiment 



data on the 5D5/2 state lifetime. 

The measurement of Ref. [ll| gave a lifetime of the 
combination of 5D^/2 and 5D5/2 states. According to 
our results, the lifetime of the whole 5d manifold should 
sit between 238.5 ns and 246.3 ns, which is consistent 
with the result of Ref. 



VI. CONCLUSION 

In this work, we have measured the lifetime of the Rb 
51)5/2 state obtaining T5D3/2 — 246.3 ± 1.6 ns, and we 
have improved the precision of the 5D^/2 lifetime result, 
T5D^/2 = 238.5 ± 2.3 ns, by a factor of 25. Our results 
have enough precision to confirm the improvement of the 
scaled all-order method [§]. 

Our experimental condition present several advantages 
over previous experimental efforts: the new laser tech- 
nique which allows us to optically excite the atoms to 
the desired states directly, the lock technique of both 
lasers in the two-step transition which makes the tran- 
sition much more efficient and the use of a MOT which 
simplifies the systematic studies. 

Our future work, the search of anapole moment in al- 
kali atoms [30l |. relies on high precision atomic structure 
calculations to extract the fundamental information from 
the measurement. We hope the advance on the experi- 
mental value of the lifetime will contribute to new im- 
provements on the precision of theoretical calculations. 
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